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INTRODUCTION 



Quantum Gravity should be the fusion between Quantum Mechanics and General Rel- 
ativity. When they are disconnected these theories work pretty good in their respective 
domains. However, when one tries to merge them in a single theory, one discovers that 
this last one is non renormalizable and therefore its predictive power is lost. Recently, 
Hofava proposed a modification of Einstein gravity motivated by the Lifshitz theory in 
solid state physics[l][2]. Such modification allows the theory to be power-counting ul- 
traviolet (UV)-renormalizable and should recover general relativity in the infrared (IR) 
limit. Nevertheless Hofava-Lifshitz (HL) theory is non-co variant, at least in its original 
formulation 1 . Indeed, in this approach space and time exhibit Lifshitz scale invariance 
of the form 

f-^andjc'-^ (1) 

with z > 1 • Z is called the dynamical critical exponent. The breaking of the 4D diffeomor- 
phism invariance allows a different treatment of the kinetic and potential terms for the 
metric: from one side the kinetic term is quadratic in time derivatives of the metric, from 
the other side the potential has high-order space derivatives. In particular the UV behav- 



See Ref.[3]. for recent progress to obtain a covariant theory. 



ior is dominated by the square of the Cotton tensor of the 3D geometry by means ofafc 6 
contribution to the propagator leading to a renormalizable power-counting theory. The 
original HL theory is based on two assumptions - detailed balance and projectability[4]. 
The projectability condition is a weak version of the invariance with respect to time 
reparametrization and therefore to the Wheeler-DeWitt (WDW) equation[5]. In this con- 
tribution, we do not want to discuss the appealing properties of the HL theory, but rather, 
we would like to investigate the original WDW equation with a look to possible defor- 
mations induced by the spacetime itself. Therefore no matter fields will be included in 
such a discussion. We have in mind two interesting cases: Non Commutative Geometry 
and Rainbow's Gravity. Both of these cases are modifications of the original space time 
structure. For this reason we expect some interesting consequences on the computation 
of some observables. In this paper, we are interested to see what are the effects of such 
distortions on Zero Point Energy (ZPE). A partial promising answer has been obtained 
in Ref.[6], where the UV divergences ogf a Schwarzschild black hole are kept under 
control by an appropriate choice of the Rainbow's metric. Units in which h = c = k = 1 
are used throughout the paper. 



THE WDW EQUATION 

The WDW equation was originally introduced by Bryce DeWitt as an attempt to quantize 
General Relativity in a Hamiltonian formulation. It is described by 



{2K)G ljkl ^K kl -^R-2A) 



m = o (2) 



and it represents the quantum version of the classical constraint which guarantees the 
invariance under time reparametrization. Gyva is the super-metric, n l i is the super- 
momentum, 3 ^ is the scalar curvature in three dimensions and A is the cosmological 
constant, while K = SnG with G the Newton's constant. In this way, the WDW equation 
is written in its most general form. The main reason to use such an equation to discuss 
renormalization problems is related to the possibility of formally re-writing the WDW 
equation as an expectation value computation 2 [7]. Indeed, if we multiply Eq.(2) by 
m* [gij\ and functionally integrate over the three spatial metric gij we find 



1 J9[g tJ ] *F [ gij ] J E dhA^ [ 8i j] _ 1 (V\kd 3 xk L \V) _ A 



(3) 



V /#fey]¥*[*y]<r[*y] " V (^> ~ K 
In Eq.(3) we have also integrated over the hypersurface E and we have defined 

V = J^d 3 x^g- (4) 

as the volume of the hypersurface E with 

A Z = (2k) G ijkl ^7t kl - y/?R/ (2k) . (5) 



See also Ref.[8] for an application of the method to a f(R) theory. 



In this form, Eq.(3) can be used to compute ZPE provided that A/jc be considered as 
an eigenvalue of A^. In particular, Eq.(3) represents the Sturm-Liouville problem as- 
sociated with the cosmological constant. To solve Eq.(3) is a quite impossible task. 
Therefore, we are oriented to use a variational approach with trial wave functionals. 
The related boundary conditions are dictated by the choice of the trial wave functionals 
which, in our case are of the Gaussian type. Different types of wave functionals cor- 
respond to different boundary conditions. The choice of a Gaussian wave functional is 
justified by the fact that ZPE should be described by a good candidate of the "vacuum 
state". To fix ideas, we choose the following form of the metric 

dr 2 

ds 2 = —N 2 (r) dt 2 + ^-j + r 2 (dG 2 + sin 2 6 d(j) 2 ) , (6) 

where b (r) is subject to the only condition b(r t ) = r t . As a first step, we decompose the 
gravitational perturbation in such a way to obtain the graviton contribution enclosed in 
Eq.(3). 



Extracting the graviton contribution 

We can gain more information if we consider gjj = gij + hij, where gjj is the back- 
ground metric and hij is a quantum fluctuation around the background. Thus Eq.(3) can 
be expanded in terms of hij. Since the kinetic part of Ax is quadratic in the momenta, 
we only need to expand the three- scalar curvature / d 3 x^/g 3 R up to the quadratic order. 
However, to proceed with the computation, we also need an orthogonal decomposition 
on the tangent space of 3-metric deformations [9, 10]: 

h,j = ±(o + 2V.Z)g tJ +{l4) iJ + htj. (7) 
The operator L maps into symmetric tracefree tensors 

(L^).. = V^ + V^-^ 7 (V^), (8) 

hfj is the traceless-transverse component of the perturbation (TT), namely g'^hfj = 0, 
V'hjj = and h is the trace of hij. It is immediate to recognize that the trace element 
a = h — 2 (V • £,) is gauge invariant. If we perform the same decomposition also on the 
momentum 7t lJ ', up to second order Eq.(3) becomes 







A^ + A| + A£ 


(2) 
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Concerning the measure appearing in Eq. (3), we have to note that the decomposition (7) 
induces the following transformation on the functional measure f^/z,-j — > ^hf-S!t,iS>aJ\, 



where the Jacobian related to the gauge vector variable is 



det 



Ag lj + ^V'Vj-R'j 



(10) 



This is nothing but the famous Faddev-Popov determinant. It becomes more transparent 
if % a is further decomposed into a transverse part £j with V a <^J = and a longitudinal 

part t}a with £,1 = V fl \]/, then J\ can be expressed by an upper triangular matrix for certain 
backgrounds (e.g. Schwarzschild in three dimensions). It is immediate to recognize that 
for an Einstein space in any dimension, cross terms vanish and 7i can be expressed by a 
block diagonal matrix. Since detA5 = detAdet.8, the functional measure @hij factorizes 
into 

i 



®hij= (det Ay) 1 ( det 



@hjj9£ T &\i/ 



(11) 



with y^v) = ^8 lJ ~R' J ac ti n g on transverse vectors, which is the Faddeev-Popov de- 
terminant. In writing the functional measure @hij, we have here ignored the appearance 
of a multiplicative anomaly[l 1]. Thus the inner product can be written as 

j @pV>* [hj^V* [£ r ]¥*[a]¥[/^]¥[£ r ]¥[(7], (12) 



where 



@p = mij$£, T @o (det Ay) 5 det 



-A 2 +V^V ; - 



(13) 



Nevertheless, since there is no interaction between ghost fields and the other compo- 
nents of the perturbation at this level of approximation, the Jacobian appearing in the 
numerator and in the denominator simplify. The reason can be found in terms of con- 
nected and disconnected terms. The disconnected terms appear in the Faddeev-Popov 
determinant and these ones are not linked by the Gaussian integration. This means that 
disconnected terms in the numerator and the same ones appearing in the denominator 
cancel out. Therefore, Eq.(9) factorizes into three pieces. The piece containing is 
the contribution of the transverse-traceless tensors (TT): essentially is the graviton con- 
tribution representing true physical degrees of freedom. Regarding the vector term A£, 
we observe that under the action of infinitesimal diffeomorphism generated by a vector 
field £ ; , the components of (7) transform as follows[9] 



h 



h + 2V-^ hj 



h 



(14) 



The Killing vectors satisfying the condition V£j + Vjfy = 0, do not change hy, and 
thus should be excluded from the gauge group. All other diffeomorphisms act on 
nontrivially. We need to fix the residual gauge freedom on the vector The simplest 
choice is = 0. This new gauge fixing produces the same Faddeev-Popov determinant 
connected to the Jacobian J\ and therefore will not contribute to the final value. We are 



left with 
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Note that in the expansion of J L d 3 Xy/gR to second order, a coupling term between the 
TT component and scalar one remains. However, the Gaussian integration does not allow 
such a mixing which has to be introduced with an appropriate wave functional. Extract- 
ing the TT tensor contribution from Eq. (3) approximated to second order in perturbation 
of the spatial part of the metric into a background term gij, and a perturbation h t j, we 
get 



where 



A L /r 



A L /r 



. -AR^j+lRhjj 



(16) 



(17) 



is the modified Lichnerowicz operator and A/ls the Lichnerowicz operator defined by 

(Al% = Ahij - 2R ikjl h kl +R lk h) + R jk h k A = -V a V a . (18) 

Qijkl re p resen t s the inverse DeWitt metric and all indices run from one to three. Note that 
the term —4R k h k L j+ 3 Rhjj disappears in four dimensions. The propagator K 1 - (x,x) iakl can 
be represented as 



*L T)± (?)*g )x (?) 

2X(x) 



(19) 



.(*)- 



where h^'^ (it) are the eigenfunctions of A L . x denotes a complete set of indices and 
X ( t) are a set of variational parameters to be determined by the minimization of Eq. ( 1 6) . 

I 

The expectation value of is easily obtained by inserting the form of the propagator 
into Eq.( 16) and minimizing with respect to the variational function X (x) . Thus the total 
one loop energy density for TT tensors becomes 



A 
87TG 



= ~iL V^m+vWm 



(20) 



The above expression makes sense only for cof (x) > 0, where (0 ; are the eigenvalues 
of A L . In the next section, we will explicitly evaluate Eq.(20) for a background of 
spherically symmetric type. 



ONE LOOP ENERGY DENSITY: CONVENTIONAL 
REGULARIZATION AND RENORMALIZATION 



The reference metric (6) can be cast into the following form 

ds 2 = —N 2 (r (x))dt 2 + dx 2 + r 2 (x) (d9 2 + sin 2 9 d(j) 2 ) , 



(21) 



where 

dr 

dx = ±^^^=. (22) 



Specific examples are 



b( r ) = ^ r 3 ; £(r) = -^r 3 and b(r)=2MG. (23) 

However, we would like to maintain the form of the line element (21) as general as 
possible. With the help of Regge and Wheeler representation [12], the Lichnerowicz 
operator ( Aj/r 1 ) . . can be reduced to 



'j 



d L 1(1+1) 2 , . 



fi(x) = (0}, fi(x) 1=1,2 , (24) 



where we have used reduced fields of the form f\ [x] = F\ (x) / r and where we have 
defined two r-dependent effective masses m\ (r) and m\ (r) 

m?(r)=6(i_^)) + ^( r) _^ (r ) 

(r = r(x)). (25) 

In order to use the W.K.B. method considered by 't Hooft in the brick wall problem[13], 
from Eq.(24) we can extract two r-dependent radial wave numbers 

kf(r^CO iM )=(Dl M - l -^^--mj(r) /=1,2 . (26) 

Then the counting of the number of modes with frequency less than ft), is given approx- 
imately by 

/"'max 

g(cOi)= Vi(l,aH)(2l + l)dl. (27) 
Jo 

V, (/, a>i) is the number of nodes in the mode with (/, ft),), such that (r = r(x)) 

1 /'+ 00 / 

V/ (/, ft),) = - / dxJkf (r, /, fti,) . (28) 

Here it is understood that the integration with respect to x and / max is taken over those 
values which satisfy kf (r, /, ft),) > 0. With the help of Eqs.(27, 28), Eq.(20) becomes 

A 1 3^ f+ ao ds(cOi) 



87rG K f^\JQ dC0i 



This is the one loop graviton contribution to the induced cosmological constant. The 
explicit evaluation of Eq.(29) gives 



A 
SkG 



1 2 r+°° I 

= Pl + p 2 = - _ £ cof^Cof - mj (r)dCOi, (30) 



where we have included an additional 4k coming from the angular integration. The use 
of the zeta function regularization method to compute the energy densities p\ and p 2 
leads to 



m. 



(r) 



64k 2 



1 



+ ln 



4/r 



m 



i = 1,2 



(3D 



where we have introduced the additional mass parameter \i in order to restore the correct 
dimension for the regularized quantities. Such an arbitrary mass scale emerges unavoid- 
ably in any regularization scheme. The renormalization is performed via the absorption 
of the divergent part into the re-definition of a bare classical quantity. Here we have two 
possible choices: the induced cosmological constant A or the gravitational Newton con- 
stant G. If we decide to absorb the divergence with the help of the cosmological constant 
A, we have to separate it into a bare cosmological constant Ao and a divergent quantity 
A div , where 

Gml (r) 



div 



A ' e32k 2 ' 

and the remaining finite value for the cosmological constant reads 



Ao 
%kG 



= ( Pl Gu) +p 2 (M)) = pJffM = ^ln ( 



m: 



(32) 



(33) 



THE EXAMPLE OF NON COMMUTATIVE THEORIES 

Non Commutative theories provide a powerful method to naturally regularize divergent 
integrals appearing in Eq.(30). Basically, the number of states is modified in the follow- 
ing way[16] 

d 3 xd 3 k , d 3 xd 3 k ( Q, 2 2/N A . in 

dn = 5- dni = ^-exp — - ( cof „, - mf (r) ) , i = 1,2. (34) 

(2k) 3 (2k) 3 V 4 v l > nl 1 W V 

This deformation corresponds to an effective cut off on the background geometry (21). 
The UV cut off is triggered only by higher momenta modes > 1/ y/0 which propagate 
over the background geometry. The virtue of this kind of deformation is its exponential 
damping profile, which encodes an intrinsic nonlocal character into fields fi(x) . Plugging 
(28) into (27) and taking account of (34), the number of modes with frequency less than 
COi,i= 1,2 is given by 

1 f+°° f/max HZ Kl+l) " / G o\ 

8(0*) = - J ^ dx j o (2/ + 1) y afa - - m] (r) exp j dl (35) 



and the induced cosmological constant becomes 



1 



A 



r j{^-ml{r))\-^ 2 -^)d<0 
J V m o( r) 



+ 



[ + °° sj \co 2 + m 2 (r)) 3 e-^ a2+m ^)dCO 
J 



(36) 



where an integration by parts in Eq.(29) has been done. By further developing the 
calculations we find 



A 



1 



where Kq (x) and (x) are the modified Bessel function and 

m 2 (r) 6 



x = 



(38) 



The asymptotic properties of (37) show that the one loop contribution is everywhere 
regular. Indeed, we find that when x — > +°°, 



A 



Conversely, when x — > 0, we obtain 



\^^[3+{Sx 2 + 6x + 3)cxp(-x)] -+0. (39) 



A 
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, , 7 3 fx 



3^ 2 2 



(40) 



a finite value for A. Note that expression (37) can be used when the background satisfies 
the relation 



m (r) = m l (r) = -m 2 (r) . 



(41) 



Examples of metrics satisfying relation (41) are the Schwarzschild, Schwarzschild-de 
Sitter (SdS) and Schwarzschild- Anti de Sitter (SAdS) backgrounds 3 . The other interest- 
ing cases, namely de Sitter and Anti-de Sitter are described by 



m \ ( r ) = m\ (r) = m\ (r) , 



(42) 



leading to 



A 



1 



%kG 3n 2 



[ + °° J{o} 2 -m 2 {r)fe-^ m2 - m ^)dco 
W m o( r ) 



Usually for such geometries, relation (41) is satisfied in a region close to the throat. 



1 



6n 2 \e 

The asymptotic expansion of Eq.(43) leads to 

A 1 



^(l-^(|) + ^(|))expg). (43) 



4\ 2 3 



6^ \ ?J 8Vx^°' 



when x — > oo. On the other hand, when x — > 0, one gets 

2 



A 



1 
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i.e. a finite value of the cosmological term 



3^ 2 2 ' 



(44) 



(45) 



RAINBOW'S GRAVITY AT WORK 



In recent years, there has been a proposal on how the fundamental aspects of special 
relativity can be modified at very high energies. This modification has been named 
Doubly Special Relativity (DSR)[17]. One of its effects is that the usual dispersion 
relation of a massive particle of mass m is modified into the following expression 

E 2 g\ (E/Epi) - p 2 g 2 2 (E/Epi) = m 2 , (46) 

where gi (E/Epi) and g2 (E/Epi) are two functions which have the following property 

lim gl {E/Epi) = \ and lim g 2 {E/E Pl ) = \. (47) 

Thus, the usual dispersion relation is recovered at low energies. This simple assumption 
has a deep impact also when the background is curved. For example, if we examine the 
Schwarzschild metric, the analysis of Magueijo and Smolin[18] shows that the energy- 
momentum tensor and the Einstein equations are replaced by a one parameter family of 
equations 

G>v (E) = SkG(E) T^ v (E) +g MV A(£) , (48) 

where G (E) is an energy dependent Newton's constant, defined so that G (0) is the phys- 
ical Newton's constant. Similarly we have an energy dependent cosmological constant 
A (E). A solution of the modified Einstein 's Field Equations for a metric of the form (6) 
is represented by a "rainbow metric" whose line element reads 

ds2 = _N^dt 2 _ ^Jr 2 ? 2 s . n2 2 



We expect the functions g\ (E/Epf) and g 2 (E/Epi) modify the UV behavior in the same 
way as GUP and Noncommutative geometry do, respectively. In presence of Rainbow's 



Gravity, we find that Eq.(2) 4 becomes 
\g\ (E) 



R- 



2A C 

sW) 



_&{E)~ lJtt " " 2K g2 (E) 

and, consequently Eq.(3) changes into 

g*(E)(V\f L d 3 xA z \V) _ A 
V 0F|Y) ?c' 

where 



¥ = 



(50) 



Az = (2k): 



'-G ijkl n^n kl - 



(51) 



(52) 



gl(E)~^~ " (2K) g2 (E) 

Of course, Eqs.(50,51) and (52) reduce to the ordinary Eqs.(2,3) and (5) when 
E/Epi — > 0. By repeating the procedure of subsection , we find that the TT tensor 
contribution of Eq. (5 1 ) is 



S\{E) 
4V 



J^d 3 x\J}& jkl 



(2k) 



g\{E) 



K 



( x i x )iikl + 



1 



(53) 



and the total one loop energy density becomes 
A -^81 (E) 82(E) 



E 2 (x) + Je 2 (x) 



(54) 



Since Eq.(24) is modified by the "rainbow metric", the two r-dependent radial wave 
numbers (26) become 



kf{r,l,(d ijn i) 



1(1 + 1) 



8l(E) 



— mf(r) i = 1, 2 



(55) 



and with the help of Eqs.(27, 28), Eq.(20) reduces to 
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i=l 



+°° d 
E igl (E)g 2 (E) — 
dEi 



E} 



\ \8\iE) 



- mf (r) dEi 



(56) 



where E* is the value which annihilates the argument of the root. In the previous 
equation we have assumed that the effective mass does not depend on the energy E. 
To further proceed, we choose a form of g\(E/Ep) and g 2 (E/Ep) which allows a 
comparison with the results obtained with a Noncommutative geometry computation 
expressed by Eq.(36). If we fix 



gx (E/E P )=exp(-a-^) 



and g 2 (E/E P ) = l, 



(57) 



Details of the calculation related to the whole section can be found in Ref.[19] 



with a > 0, Eq.(56) becomes 
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/ E 2 exp(-a^)J(E 2 -m 2 (r))dE 



E 2 



/■+ 00 „ E 2 l 

+ E 2 zM-a-2)\j{E 2 + m 2 (r))dE 



(58) 



(59) 



where we have used relation (41). Nevertheless we have to observe that even if the final 
result is finite and background dependent, the induced cosmological constant computed 
with the choice (57) will be always negative. When we compare this result with that 
obtained in Noncommutative geometry, namely Eq.(36), we find that the negativity is 
principally due the Rainbow's functions that do not completely enter in the counting of 
nodes like in expression (35). Therefore the pure "Gaussian" choice (57) can not give a 
positive induced cosmological constant. We are thus led to choose 



gi (E/Ep) 



E \ E 2 

1+/s _j exp( _ a _ 



and g 2 (E/E P ) = l, 



(60) 



with a > and /3 e K. Plugging parametrization (57) into expression (56), one gets 
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3x x 



3v^ 



2a 2 



, sinh (ax 2 ) + -^- T cosh (ax 2 ) + 
«2 2a? 2a? 



ax 



r - JL) e ax2 erf (Va~x) 



. (61) 

where x = \Jm^ (r) /Ep and where we have used expression (41). The asymptotic 
expansion for large x is 



A ( 2 ^« 3/2 + v^ a2 ) x _ 8/3a 5 / 2 + 3^a 3 3 16/3a 7 / 2 + 5y/jm 4 

\7iG~ Aa 1 / 2 16a 11 / 2 * + 128 a l5 / 2 x 3 



SnG Aa 1 / 2 
while for small x, one gets 
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16a 11 / 2 * 
4a 5 / 2 + 3^/3 a 2 



4a 9 / 2 



+ 0(x 3 ) 



If we set 



P = 



'an 



-0(x-*), 
(62) 

(63) 

(64) 



then the linear divergent term of the asymptotic expansion (62) disappears and Eq. (61) 
vanishes for large x, while for small x we get 

A 37T-8 



SnG Sa 2 



+ 0(x 3 ) 



(65) 



where we have used the result of expansion (63). It is possible to show that with choice 
(64), the induced cosmological constant is always positive. 



CONCLUSIONS 



In this contribution we have discussed how modifications of some basic points of Gen- 
eral Relativity lead to finite contributions of ZPE without invoking neither a regulariza- 
tion nor a renormalization procedure. In particular, we have considered the effect that 
a NCG has on the counting of states with a deep modification of the measure in phase 
space. Such a modification introduces a Gaussian damping with a natural length: the 
Noncommutative length 9. The final result is a finite ZPE, interpreted as a Cosmologi- 
cal Constant, regular at each point of the spacetime. On the other hand we have MDR's 
which provide an alternative way to keep under control UV divergences [6]. The appli- 
cation of MDR's is done in terms of a "rainbow metric" which, unfortunately, predicts a 
negative cosmological constant when one considers the same damping factor suggested 
from NCG. Nevertheless, by introducing appropriate variations of the original Gaussian 
proposal, positive contributions can be obtained. 
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